
********years of education over periods******
use graphs_years_education.dta,clear

***dummy for female****
gen d_coefficient_f_p= d_coefficient_f + 1.96* d_std_f
gen d_coefficient_f_l= d_coefficient_f - 1.96* d_std_f

sencode name,gen(cohort)

***dummy for female****
twoway (rcap d_coefficient_f_p d_coefficient_f_l cohort in 1/4, xlab(1 (1) 4, val angle(45)) ) || scatter  d_coefficient_f cohort in 1/4, graphregion(color(white)) legend(label(1 "95% CI") label(2 "true coefficients")) ytitle("treatment effects")



